Autosomal Dominant Familial Dyskinesia and Facial Myokymia Single Exome Sequencing Identifies a Mutation in Adenylyl Cyclase 5

Background: Familial dyskinesia with facial myokymia (FDFM) is an autosomal dominant disorder that is exacerbated by anxiety. In a 5-generation family of German ancestry, we previously mapped FDFM to chromosome band 3p21-3q21. The 72.5-Mb linkage region was too large for traditional positional mutation identification. Objective: To identify the gene responsible for FDFM by exome resequencing of a single affected individual. Participants: We performed whole exome sequencing in 1 affected individual and used a series of bioinformatic filters, including functional significance and presence in dbSNP or the 1000 Genomes Project, to reduce the number of candidate variants. Co-segregation analysis was performed in 15 additional individuals in 3 generations. Main Outcome Measures: Unique DNA variants in the linkage region that co-segregate with FDFM. Results: The exome contained 23 428 single-nucleotide variants, of which 9391 were missense, nonsense, or splice site alterations. The critical region contained 323 variants, 5 of which were not present in 1 of the sequence databases. Adenylyl cyclase 5 (ADCY5) was the only gene in which the variant (c.2176G>A) was co-transmitted perfectly with disease status and was not present in 3510 control white exomes. This residue is highly conserved, and the change is nonconservative and predicted to be damaging. Conclusions: ADCY5 is highly expressed in striatum. Mice deficient in Adcy5 develop a movement disorder that is worsened by stress. We conclude that FDFM likely results from a missense mutation in ADCY5. This study demonstrates the power of a single exome sequence combined with linkage information to identify causative genes for rare autosomal dominant mendelian diseases

Copy Number Variation Analysis in Single-Suture Craniosynostosis: Multiple Rare Variants Including RUNX2 Duplication in Two Cousins With Metopic Craniosynostosis

Little is known about genes that underlie isolated single-suture craniosynostosis. In this study, we hypothesize that rare copy number variants (CNV) in patients with isolated single-suture craniosynostosis contain genes important for cranial development. Using whole genome array comparative genomic hybridization (CGH), we evaluated DNA from 186 individuals with single-suture craniosynostosis for submicroscopic deletions and duplications. We identified a 1.1 Mb duplication encompassing RUNX2 in two affected cousins with metopic synostosis and hypodontia. Given that RUNX2 is required as a master switch for osteoblast differentiation and interacts with TWIST I, mutations in which also cause craniosynostosis, we conclude that the duplication in this family is pathogenic, albeit with reduced penetrance. In addition, we find that a total of 7.5% of individuals with single-suture synostosis in our series have at least one rare deletion or duplication that contains genes and that has not been previously reported in unaffected individuals. The genes within and disrupted by CNVs in this cohort are potential novel candidate genes for craniosynostosis. (C) 2010 Wiley-Liss, Inc

Personalized copy number and segmental duplication maps using next-generation sequencing

Despite their importance in gene innovation and phenotypic variation, duplicated regions have remained largely intractable owing to difficulties in accurately resolving their structure, copy number and sequence content. We present an algorithm (mrFAST) to comprehensively map next-generation sequence reads, which allows for the prediction of absolute copy-number variation of duplicated segments and genes. We examine three human genomes and experimentally validate genome-wide copy number differences. We estimate that, on average, 73-87 genes vary in copy number between any two individuals and find that these genic differences overwhelmingly correspond to segmental duplications (odds ratio = 135; P < 2.2 x 10(-16)). Our method can distinguish between different copies of highly identical genes, providing a more accurate assessment of gene content and insight into functional constraint without the limitations of array-based technology

Using the posterior distribution of deviance to measure evidence of association for rare susceptibility variants

Aitkin recently proposed an integrated Bayesian/likelihood approach that he claims is general and simple. We have applied this method, which does not rely on informative prior probabilities or large-sample results, to investigate the evidence of association between disease and the 16 variants in the KDR gene provided by Genetic Analysis Workshop 17. Based on the likelihood of logistic regression models and considering noninformative uniform prior probabilities on the coefficients of the explanatory variables, we used a random walk Metropolis algorithm to simulate the distributions of deviance and deviance difference. The distribution of probability values and the distribution of the proportions of positive deviance differences showed different locations, but the direction of the shift depended on the genetic factor. For the variant with the highest minor allele frequency and for any rare variant, standard logistic regression showed a higher power than the novel approach. For the two variants with the strongest effects on Q1 under a type I error rate of 1%, the integrated approach showed a higher power than standard logistic regression. The advantages and limitations of the integrated Bayesian/likelihood approach should be investigated using additional regions and considering alternative regression models and collapsing methods

The Facioscapulohumeral muscular dystrophy region on 4qter and the homologous locus on 10qter evolved independently under different evolutionary pressure

BACKGROUND: The homologous 4q and 10q subtelomeric regions include two distinctive polymorphic arrays of 3.3 kb repeats, named D4Z4. An additional BlnI restriction site on the 10q-type sequence allows to distinguish the chromosomal origin of the repeats. Reduction in the number of D4Z4 repeats below a threshold of 10 at the 4q locus is tightly linked to Facioscapulohumeral Muscular Dystrophy (FSHD), while similar contractions at 10q locus, are not pathogenic. Sequence variations due to the presence of BlnI-sensitive repeats (10q-type) on chromosome 4 or viceversa of BlnI-resistant repeats (4q-type) on chromosome 10 are observed in both alleles. RESULTS: We analysed DNA samples from 116 healthy subiects and 114 FSHD patients and determined the size distributions of polymorphic 4q and 10q alleles, the frequency and the D4Z4 repeat assortment of variant alleles, and finally the telomeric sequences both in standard and variant alleles. We observed the same frequency and types of variant alleles in FSHD patients and controls, but we found marked differences between the repeat arrays of the 4q and 10q chromosomes. In particular we detected 10q alleles completely replaced by the 4q subtelomeric region, consisting in the whole set of 4q-type repeats and the distal telomeric markers. However the reciprocal event, 10q-type subtelomeric region on chromosome 4, was never observed. At 4q locus we always identified hybrid alleles containing a mixture of 4q and 10q-type repeats. CONCLUSION: The different size distribution and different structure of 10q variant alleles as compared with 4q suggests that these loci evolved in a different manner, since the 4q locus is linked to FSHD, while no inheritable disease is associated with mutations in 10qter genomic region. Hybrid alleles on chromosome 4 always retain a minimum number of 4q type repeats, as they are probably essential for maintaining the structural and functional properties of this subtelomeric region. In addition we found: i) several instances of variant alleles that could be misinterpreted and interfere with a correct diagnosis of FSHD; ii) the presence of borderline alleles in the range of 30–40 kb that carried a qA type telomere and were not associated with the disease

The Impact of Errors in Copy Number Variation Detection Algorithms on Association Results

The inaccuracy of copy number variation (CNV) detection on single nucleotide polymorphism (SNP) arrays has recently been brought to attention. Such high error rates will undoubtedly have ramifications on downstream association testing. We examined this effect for a wide range of scenarios and found a noticeable decrease in power for error rates typical of CNV calling algorithms. We compared power using CNV calls to the log relative ratio and found the latter to be superior when error rates are moderate to large or when the CNV size is small. It is our recommendation that CNV researchers use intensity measurements as an alternative to CNV calls in these scenarios

Family- and population-based designs identify different rare causal variants

Both family- and population-based samples are used to identify genetic variants associated with phenotypes. Each strategy has demonstrated advantages, but their ability to identify rare variants and genes containing rare variants is unclear. To compare these two study designs in the identification of rare causal variants, we applied various methods to the population- and family-based data simulated by the Genetic Analysis Workshop 17 with knowledge of the simulated model. Our results suggest that different variants can be identified by different study designs. Family-based and population-based study designs can be complementary in the identification of rare causal variants and should be considered in future studies

Single-Tissue and Cross-Tissue Heritability of Gene Expression Via Identity-by-Descent in Related or Unrelated Individuals

Family studies of individual tissues have shown that gene expression traits are genetically heritable. Here, we investigate cis and trans components of heritability both within and across tissues by applying variance-components methods to 722 Icelanders from family cohorts, using identity-by-descent (IBD) estimates from long-range phased genome-wide SNP data and gene expression measurements for ∼19,000 genes in blood and adipose tissue. We estimate the proportion of gene expression heritability attributable to cis regulation as 37% in blood and 24% in adipose tissue. Our results indicate that the correlation in gene expression measurements across these tissues is primarily due to heritability at cis loci, whereas there is little sharing of trans regulation across tissues. One implication of this finding is that heritability in tissues composed of heterogeneous cell types is expected to be more dominated by cis regulation than in tissues composed of more homogeneous cell types, consistent with our blood versus adipose results as well as results of previous studies in lymphoblastoid cell lines. Finally, we obtained similar estimates of the cis components of heritability using IBD between unrelated individuals, indicating that transgenerational epigenetic inheritance does not contribute substantially to the “missing heritability” of gene expression in these tissue types

The Role of AGG Interruptions in the Transcription of FMR1 Premutation Alleles

Fragile X associated disorders are caused by a premutation allele in the fragile X mental retardation 1 gene (FMR1) and are hypothesized to result from the toxic effect of elevated levels of expanded FMR1 transcripts. Increased levels of FMR1 mRNA have indeed been reported in premutation carriers; however the mechanism by which expanded alleles lead to elevated levels of FMR1 mRNA in premutation carriers is unknown. Within the CGG repeat tract AGG interruptions are found, generally 1–3 present in normal/intermediate alleles (6–54 CGG repeats) and usually 0–1 in premutation alleles (55–200 CGG repeats). They are present at specific locations, generally occurring after 9 or 10 uninterrupted CGG repeats [(CGG)9AGG(CGG)9AGG(CGG)n]. We evaluated both the number of AGG interruptions and the resulting length of the uninterrupted 3′ CGG repeat pure tract in premutation alleles derived from two large cohorts of male and female carriers to determine whether the presence of AGG interruptions or the length of a pure stretch of CGG repeats influence the levels of FMR1 mRNA in blood. Our findings indicate that neither the number of AGG interruptions, nor their position along the CGG tract have a significant affect on mRNA levels in premutation carriers. We also, as expected based on previous findings, observed a highly significant correlation between CGG repeat number (as both total length and length of pure CGG stretch) and FMR1 mRNA expression levels, in both males and females. Importantly, we did not observe any significant difference in FMR1 mRNA levels in premutation carriers based on age

An evaluation of different meta-analysis approaches in the presence of allelic heterogeneity

Meta-analysis has proven a useful tool in genetic association studies. Allelic heterogeneity can arise from ethnic background differences across populations being meta-analyzed (for example, in search of common frequency variants through genome-wide association studies), and through the presence of multiple low frequency and rare associated variants in the same functional unit of interest (for example, within a gene or a regulatory region). The latter challenge will be increasingly relevant in whole-genome and whole-exome sequencing studies investigating association with complex traits. Here, we evaluate the performance of different approaches to meta-analysis in the presence of allelic heterogeneity. We simulate allelic heterogeneity scenarios in three populations and examine the performance of current approaches to the analysis of these data. We show that current approaches can detect only a small fraction of common frequency causal variants. We also find that for low-frequency variants with large effects (odds ratios 2–3), single-point tests have high power, but also high false-positive rates. P-value based meta-analysis of summary results from allele-matching locus-wide tests outperforms collapsing approaches. We conclude that current strategies for the combination of genetic association data in the presence of allelic heterogeneity are insufficiently powered